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1 Introduction 

We will consider the system of linear congruenses, 

auxi + ai2X2 + • • ■ + ainXn = bi (mod m) 
a2ixi + 0222:2 + • ■ • + a2nXn = &2 (mod m) 

(1) 

^ttnixi + an2X2 H h a„„2:„ = bn (mod m), 

where to, n are positive integers, and all a^j, bi are integers. We are interested 
in finding an expression for the number of solutions to this system. Let A = 
(ttij) denote the coefficient matrix. It is well know that this system has a 
unique solution if and only if det{A) and m are relatively prime, see for example 
[8j. But, how many solutions do we have if det{A) and m have a common 
divisor greater than 1? The authors became interested in this question when 
they where working on multidimensional p-adic monomial dynamical systems 
[7]. The problem was solved by Butson and Stewart, [I], by rewriting the 
system into Smith normal form. The solvability of linear congruence systems and 
algorithms for finding a solutions have been of interest for many mathematicians 
and computer scientists over the years. Example of other contributions to the 
problem on solving systems of linear congruences can be found in [5], [5], [3], [1] 
and [g. 

Compared to Butson and Stewart we will use a more direct method. We 
use Gaussian elimination with successive reduction and the Chinese remainder 
theorem instead of the Smith normal form. We will find a different formula for 
the number of solutions than in [T]. The algorithm used in the proof of the 
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formula can with small changes also be used to find all incongruent solutions to 
the system, since we only have used elementary methods. 

The paper is organized as follows: In Section [5] we give definitions and no- 
tations together with some theorems about solvability of systems of linear con- 
gruences. Section [3] is the main section in which we derive the formula for the 
number of solutions to homogeneous systems. Applications to inhomogeneous 
systems are mentioned in Section In Section [S] we present an algorithm in 
pseudo code for calculating the number of solutions. The algorithm also gives 
us all the solutions. In Section [S] we discuss generalizations of our results. 

2 Notations 

Two solutions of the linear congruence system ([T]) are said to be incongruent 
modulo m if they differ at least in one coordinate modulo m. We want to 
find the number of incongruent solutions modulo m the system have. Let A 
denote the coefficient matrix, x the vector of the indeterminates, and b the 
vector of the elements on the right hand side in the system, that is, A = {a,ij), 
X = (xi, X2, ■ ■ ■ ,Xn) and b = 62, . . . , &n). Then the system ([Ij can be written 
in matrix form as 

Ax = b (mod to). (2) 

Let r](A, 6, to) denote the number of incongruent solutions modulo to to the 
congruence 

Let adj(/4) denote the adjoint matrix of a square matrix A over order n. It 
is known from linear algebra that adj(/4) has the following properties: 

det(adj(/A)) = det(/A)""i 

and 

A adi{A) = adj(>^\) A = dct{A) /, 
where / is the identity matrix. 

Theorem 2.1. Let A, b and to be as above. Then 

v{A,b,m) < (det(/\),TO)". 
If {det{A),m) = 1, then ri{A,b,m) = 1. 

Proof. We multiply the congruence system Ax = b (mod m) with adj(/4) and 

get 

adj {A) Ax = adj {A) b (mod to) 

det{A) X = a,di{A) b (mod to), (3) 

which is solvable if and only if det{A) divides all elements in the vector adj(/4)6. 
If that is the case then this system has (det(/4),TO)" different solutions. If that 
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is not the case, then the original system has no solutions, since any solutions 
is also solutions to the rewritten system ([3]). When we multiply with adj(/4) it 
might happen that we introduce new solutions. This proves that ri{A, b, m) < 
(det(/\),m)". 

Assume that (det(/4), to) = 1. From ^ it follows that 

X = det(/4)"^ adj(/4) h (mod m) 

is a solution to the original system. Hence r\{A^ b, m) — 1. □ 

Theorem 2.2. Let A, b and m be as above. Assume that to, = TOi ■ • • to^, where 
the integers toi, . . . , toj, are pairwise relatively prime. Then 

r]{A, b, m) = t]{A, b, mi) • • • ri{A, b, m^). 

Hence, rj is multiplicative with respect to to. 

Proof. If a = 6 (mod to) and n divide to, then a = b (mod n). Hence any 
solution to Ax = b (mod to) is also a solution to every Ax = b (mod rui), 
where i = 1,2, ... ,k. Assume that Xi — {xn, Xi2, . . . , Xin) is a solution to the 
congruence system Ax = b (mod to^), for i = 1,2, . . . , k. Then for any fc-tuple 
{xi,X2, . . . ,Xf^) of solutions of the systems modulo to^ we construct a solution 
X = (xi, X2, . . . , x„) modulo to, by defining Xi to be the unique solution, due to 
the Chinese Remainder Theorem, modulo to, to the system 

X = xii (mod TOi) 
X = X2i (mod TO2) 

^x = xut (mod TOfc). 

Since there are r]{A,b,mi) ■ ■ ■ r}{A,b,mk) possible such fc-tuples of solutions 
modulo TOi, the theorem follows (two different fc-tuples can not generate the 
same solution modulo to). □ 

3 Homogeneous Systems 

We have already stated that we only have to consider congruence systems mod- 
ulo a prime power. So, from now on to is going to be equal to p^ , where p is a 
prime number and fc a positive integer. Let d — {A,p^), the greatest common 
divisor of all elements in the matrix A and the integer p^ , that is, 

{A,p'') = (oii,oi2,...,a„„,p''). 

Hence, d ~ p'^ for some non-negative integer e < k. 

Lemma 3.1. Let d = {A,p'') and let R = A/d. Then 

r;(A0,/) = d",7(R,0,//d). 

Note that r]{R,0, 1) = 1. 
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Proof. Let x' be a solution to Rx = (modp'^/d). Note that x' is also a 
solution to Ax = (mod p''). This is also true for all elements on the form 



where u. G Z^. This proves the theorem since jZ^I = d". □ 

Let — /d and let R be as in Lemma [3.11 We know that — 1 and 

therefore we can find some element in R that is relatively prime to p, say rn after 
a possible rearrangement of rows and columns. Note that rn is then invertible 
modulo p' . By Gaussian elimination we get the equivalent system 



riixi 



(1) 



(1) 



ri3X3 

(1) 



(1) 

''»-l,22;3 



^InXn 
' l,n-l-^ 



.(1) 



' n-l,n-l-^n 



(mod p') 
(mod p') 

(modp'), 



(4) 



where 



for i,j = 
order n 



(mod p') 



1,2,. 



, n — 1. Let /?(^^ = ('^ij^)- Hence, /?(^) is a square matrix of 



1. Since rn is invertible modulo the two systems 

Rx = (modp') and R^^^x^^'>=0 (mod p') 

has the same number of solutions. We have proved the first part of the following 
lemma. 



Lemma 3.2. With the notation above, 



have 



r,{R,0,p')^rj{R(^\0,p'). 
Moreover, we have that (det(/?),p') — (det(/?(^'),p'). 
Proof. From the properties of the determinant it follows that 
det{R) = rii(-l)i+Met(/?(i)) (mod p'). 
Since (rii,p) = 1 we have (dct(/?),p') = (det(/?(i)),p'). 
Example 3.1. Let p be a prime number, I a positive integer, and 



(5) 



□ 



R 



ni ri2 

f21 r22 



an integer matrix. Assume that {R,p) — (?'ii, fi2, ^21, ''22,p) — 1- Consider the 
linear congruence system 



riiXi + ri2X2 = (mod p') 
r2iXi + r22X2 = (mod p'). 



(6) 
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Since {R,p) = 1 one of the matrix entries must be relatively prime to p, assume 
that it is rii. Let denote the multiplicative inverse of rn modulo pK From 
the first congruence in the system ^ we have 

xi = ~r^iri2X2 (mod p'). 

Hence, xi is uniquely determined by X2 modulo . By putting this expression 
for xi in the second congruence in the system ^ we get 

{r22 - r^iri2r2i)x2 = (mod p'), 

or equivalent 

(?'ii?'22 - ?'i2?'2i)a;2 = (modp'). 

The coefficient for X2 is det(/?) and from the theory of linear congruences we 
know that this equation has (det(/?),p') incongruent solutions modulo pK We 
conclude that ?](/?, 0,p') = (det(/?),p'). 

Theorem 3.3. Let p be a prime number, I a positive integer, and A an integer 
square matrix of order n. Then 

77(AO,/)=p"'=-^ (7) 

where j — lo + li + ■ ■ ■ + In-i- The numbers k, for i = 0,1, . . . ,n — l, are defined 
recursively in the following way: 

p''=4-^, and = (/?«,p''-0, 
di 

with do = {A,p''), p^° = P^ /do and where Z?*^'' is the matrix of order n — i given 
analogously as in Q from R^^~^\ 

Proof. From Lemma l3. II we have that 

7?(Ao,/) = rfo^K^^°\o,p'''), 

where do = {^tP^) and p'" = p^ /do. If we now construct R'^'^^ from A/do in the 
same way as in the proof of Lemma 13.21 we get that 

7j{A,Q,p^) = dlin{R^^\0,p'-). 

We now do the same for the system R^^^a; = (mod p^°). If Z?^^-* = (mod p^"), 
then ri{R''^\0,p^°) = p'o("-i) and in this case 7 = Iq- Otherwise, we construct 
the matrix Z?'-^^ of order n — 2 from /?(^) in the same way as we constructed /?(^) 
before. Let di = {R^^\p'-^) and p'^ —p'-°/di. By using Lemma [3. II and Lemma 
13.21 again we get 

r/(AO,/) ^d'^d1-^r]{R'^^\0,p^'). 

If = (mod p'l) then r/(/?(2), 0,p'i) = p'l and j ^ lo + h- If we continue 
in this way we get 

77(4, 0,/) = d^dr' ■ ■ ■ dl_2v{R'^--^\0,p'-'^). 
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The matrix is of order 1, that is, a single integer r^l ^\ Let 

= (/?("-!), p"-2) = (r(r'\p'"-^) and p''^-^ = p^'^-^dn-i. 
The the congruence R^"'~^^x = (mod is equivalent to 

r[i ^'a; = (modp'""^), 
which have dn-i incongruent solutions modulo p'"-^. Hence 

7?(/?("-l),0,p''-) = rf„_i. 

We have showed that 

Let l-i = k. From dt = p''- Vp'* = it follows that 

d^d^-' ■ ■ ■ di_^dn-i = n = n p~'* = f"'"''. 

i=0 1=0 

since 7 = + /i H h L-i- This proves that r]{A, 0,p'') = □ 

Theorem 3.4. With the same notations as above, (det(y4),p"''') =p'^''~'^. 

Proof. Let =A,do = d, = p'^-i), where /?W is constructed 

as above. Let also 5^*^ = R'^^^/di. Hence is the submatrix we get from 

Gaussian elimination in S^'^ by excluding the first row and column. It is clear 
that 

dct{A) = det(/?(°)) = det(S(")). 
After Gaussian elimination and expansion along the first column we have 

det(S(°)) = det(/?(i)) (mod p^") 

and 

(det(S(°)),p'«) = (det(/?(^)),p'«). 

Hence, 

(det(4),djjp'") = (det(S(°))dJJ,dJJp'«) = (det(i?(^))d^,d^p'''). 
Observe that dQp'"> = p"'=-("-i)'o, general we have 

det(/?(^)) = rf^-^det(S(^)) 

and 

(det(S(^)),/') = (det(/?(^+i)),p'0> 
for J = 0. 1, .... ri — 1. By induction 

idem, d^... d]Z(+'p'^) = (det(/?(^"))dS • • • d]Z(+\d^ ■ ■ ■ d-Zl+^p'^), 
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for all j. For j = n ^ 1 we get 

(det(/\), d^ • • d2_2p''-) = (dct(i?("-i))d^ • • • dl_^,d^ ■ ■ ■ dl^,p'-^) 

Since p"'^ > rfg • ■ • d^_2p'""^ > p"'"'"'*' we conclude 

(det(/\),p"'=) □ 

Corollary 3.5. The number of solutions of Ax = (mod p'') is 

77(AO,/) = (det(/\),p"'^-). 

Theorem 3.6. Let A be a square integer matrix of order n, and m a positive 
integer. Assume that m — Jlfc=iPi^S where pi are different primes. Then 

N 

7j{A,0,m) = m"Y[pp' = (det(/l),m") 

i=l 

where ji is defined as in Theorem \3.3[ 

Proof. This is a direct consequence of Theorem 12.21 and Theorem 13.31 □ 
From Theorem 12.11 and Theorem 13.61 we have 

r]{A,0,m) = (det(4),m") < (det(/\), m)" 
and with equality when (det(^),m)" | det(A). 

4 Inhomogeneous Systems 

Theorem 4.1. Let p be a prime number, k a positive integer, A be a square 
integer matrix of order n, and b an integer vector of length n. The inhomogeous 
linear system 

Ax = b (mod p*^) (8) 

is solvable if and only if do \ b and dj \ b*--'' for all < j < n — 1, where the 
integers dj are given in similar way as in Theorem \3.3\ and the vectors b^-'^ is 
the last n — j elements in the right hand side of the system after the jth .step of 
the Gaussian elimination. Moreover, if the system is solvable then 

ry(Ab,/) = ?y(AO,p'=) =p"'=-^ = (det(^),p"'=). 

Proof. Let do = {A,p'^). The system ([5]) is solvable when do divide each com- 
ponent of b. Let = A/do and b*-"-* — b/do. Reduce the system in the same 
way as in Lemma 13.21 the first step in the Gaussian elimination. We get the 
inhomogeous system 

/?(i)a;(i) = (mod 
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with 71 — 1 unknowns and where 

6f)=fe|°)-5j°V,|°)(r(°))-i (modp'°) for i = 2,3,...,n. 

In order for this system to have solutions c?i = {R^^\p'-") must divide all the 
components of b*'^' . If we continue in this way we get that is solvable if and 
only if do I b and dj \ b^-'^ for all < j < n — 1. The number of solutions, if they 
exists, are the same as in homogeneous case — the backward substitution result 
in the same number of solutions in each step. □ 



5 Algorithm 

Let p be a prime number, n and k positive integers, A — (ay)„xn an integer 
matrix, and b — (bi)nxi a vector with integer entries. The following algorithm 
decides if the congruence 

Ax = b (modp'') 

is solvable, and in that case find the number of solutions, denoted rj, and the 
set X of all solutions. In the algorithm we will work with matrices on the block 
matrix form 




t = 1,2,. 1, 



where L/*^*^ is an upper triangular square matrix of order i, 6^*^ is a t x n — t 
matrix, is the zero matrix of type in — t) x i, and Z?*^*^ is a square matrix 
of order n~t. Let = A'-°\ Hence, is an upper triangular matrix. 

Let rf-' denote the element on the ith row and the jth column in Z?*^*^. Then 



Let bf' denote the ith element in the vector b*-*-*. Set 

that is, the n — t last elements in b*-*'. Hence, r'") = b^'^K Note that when we 
change /?^*' or r^*) we also change /A*^*) and b^*-*, respectively. Let U{q, n) = 
1. Initializaion 

1.1. l-i ^ k 

1.2. a <— id, the permutation which represent the interchanges of columns 

1.3. do ^ (/4,p'^), note that do = p'^" for some integer eo 

1.4. lo ^ k ~ eo, that is p'° — /d^ 

1.5. 7 ?o 

1.6. /iW ^ A/da 
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1.7. If rfo I b, then ' b/do, else goto step 5 

1.8. If da =p^ , then X ^ U{p^,n) and goto step 4, note that 7 = 

2. Gaussian elimination 
For f = 1,2, ...,n- 1 do 

2.1. Find u and v such that (r'ujr^\p) = 1 

2.2. Interchange tth and [t + u— l)th row, and tth and if + v — l)th column 
of Let denote the new matrix. 

2.3. Interchange tth and {t + u— l)th elements of b^*'~^\ Let 6^*^ denote the 
new vector. 

2.4. Swap: (7{t) ^ a{t + v-l) and (7(t + - 1) -S- <T(t) 

2.5. For i,j = t + l,t + 2, . . . ,n do 

2.5.1. agVag)-(ai:))-^(fag)modp'*- 

2.5.2. 6fV6f -(a(:V'«i^^*^^lodp'*- 
2.5.3. a^^*^ ^ 

2.6. <— {R'^*\p''*-^), note that dt = p'^* for some integer 

2.7. ^ /t_i - et, that is p'* = p'-^'^/dt 

2.8. 7 <«-7 + Zt 

2.9. /?W ^ 

2.10. If dt I r(*), then r^*) ^ r^'V^^t. else goto step 5 

2.11. If dt = _p''-S then X -i- U{p'-*-\n- t) and goto step 3 

2.12. If f = n- 1, then 

2.12.1. Xn (a^'Jr^'') modp'"-i 

2.12.2. X <^Xn +p'"-'U{d„-i, 1) modp'"-^ 

3. Backward substitution 
For s = 1,...,1 do 

3.1 y^0 

3.2 For each x = {xs+i,Xs+2, ■ ■ ■ ,Xn) G X do 

3.2.1 xs ^ (ai*i)"'[6i*) - ^ aj^xj") modp'^-^ 

3.2.2 X [xg^Xs^i, . . . that is, prepend Xg to x 

3.2.3 y ^ F U (a; + p'"--'U{ds-i,n -s + 1) mod p''-^) 

3.3 X ^ r 

4. Solvable 

4.1. For each x = {xi,X2, . . . ,a;„) G X do a; <(- , a;<^(2) , • • ■,Xa{n)) 

4.2. -s- 

4.3. Goto step 6 
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5. Not solvable 

5.1. T]i~0 

5.2. X ^ 

6. Output and finish 
6.1. Return ry and X 

Remark. After step 1.7: we have transformed 



Ax = b (mod/) to A^°^x = b^°^ (mod/). 



After step 2.10: we have transformed 

/?(t)aj = rW (mod /-I) to 



/?Wa; = rW (mod/) 



since p'* = p''-^/dt. During the backward substitution solve first modulo p'* 
and thereafter lift to modulo p''-^. 

Remark. If we arc only interested in finding the number of solutions its prefer- 
able to use the formula with the determinant. 

Example 5.1. Study the linear congruence system 



' 123a;i + 152a;2 + 28x3 + 22x4 + 144x5 = 193 
38x1 + 189x2 + 127x3 + 171x4 + 141x5 = 2 
132x1 + 232x2 + 215x3 + 22x4 = 96 

155x1 + 30x2 + 178x3 + 142x4 + 127x5 = 198 
194x1 + 171x2 + 16x3 + 24x4 + 98x5 = 162 



where 3 = 243. Any solutions are elements in the set 



(mod 3^), 



X CU = Z243 X Z243 X Z243 X Z243 X Z. 



Let 



/123 152 
38 189 



A= ^ 

155 30 
\194 171 

Set l-i = 5. We have that eo 



28 22 
127 171 



132 232 215 22 



178 142 
16 24 



id: 



0, do 

1 
1 



144\ 
141 


127 
98 / 

1, lo 



^243- 

/193\ 
2 

and b = 96 
198 
Vl62/ 

5, and 7 = Zo = 5. Initial we let 



2 3 4 5 
2 3 4 5 



Since do = 0, we have that = A and 6^°) = b. Note that R^°'> = A^^l We 
interchange row 1 and 2. After the first step in the Gaussian elimination we get 



^(1) 



/38 189 127 171 14l\ 

71 7 76 180 

151 68 157 9 

84 114 52 121 

\ 63 135 123 56 / 



and 6^^) = 



/ 2 \ 
97 
153 
241 

V139/ 
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The system given by the matrix A^^^ and vector b'^^^ have the same solutions as 
original system. Further, the new system have as many solutions as 



= ^(1) (mod 3^), 



where 



(n 7 76 180^ 

151 68 157 9 

84 114 52 121 

\ 63 135 123 56 / 



and r(i) 



/97\ 
153 
241 

Vi39y 

lo + h 



For this system we find that ei = 0, di = 1, li = 5, and ^ = lo + h = 10. One 
more step with the Gaussian elimination results in 



4(2 



Next we have to interchange column 3 and 4 since none of 36, 27 or 9 are relative 
prime to 3. Hence, the column permutations is so far given by 



1 2 3 4 5 
1 2 4 3 5 



/38 


189 


127 


171 


141\ 




( 2 \ 





71 


7 


76 


180 


and 6^2^ = 


97 








36 


122 


54 


22 








27 


10 


175 




181 


Vo 





9 


213 


218^ 







We get that 62 = Q, d2 = 1, h = 5, and 7 
the Gaussian elimination gives us that 



^(3) 



^0 + ^1+^2 = 15. The next step in 



/38 


189 


171 


127 


141\ 




( 2 \ 





71 


76 


7 


180 


and b^^^ = 


97 








122 


36 


54 


22 











36 


67 




227 


vo 








225 


56 J 




i,i99y 



Again, we have to interchange columns since both 36 and 225 are not relative 
prime to 3. That gives us the column permutation 



a = 



1 2 3 4 5 
1 2 4 3 5 



After the last step in the Gaussian elimination we have that 



/38 


189 


171 


141 


127\ 




( 2 \ 





71 


76 


180 


7 




97 








122 


54 


36 


and 


22 











67 


36 




227 


^0 











126^ 
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and 63 — 0, = 1, = 5, and j = Iq + h + h + h = 20. Finally, we get that 
64 = 2, d4 = 9, U = 3, and 7 = + ^1 + '2 + ^3 + ^ = 23. Since ^4 = 2^ = 9 
divide both 126 and 216, the system is solvable, and we have that 



Now, we do backward substitution on the linear congruence system 
r 382/1 



/38 


189 


171 


141 


127\ 




f 2 \ 





71 


76 


180 


7 


and b^^) = 


97 








122 


54 


36 


22 











67 


36 




227 


Vo 











14 yl 




I 24/ 



189y2 

71J/2 



171y3 
762/3 
122J/3 



1412/4 + 127J/5 
I8O2/4 + 72/5 
547/4 + 36j/5 
671/4 + 36j/5 
14y5 



2 
97 
22 
227 

24 



(mod 3^) 
(mod 3^) 
(mod 3^) 
(mod 3^) 
(mod 3^), 



where 



38-1 
71-1 

122-1 
67-1 
14-1 



32 
89 
2 

214 
2 



(mod 3^) 
(mod 35 ) 
(mod 3^^) 
(mod 35 ) 
(mod 33) 



and 



xi = ya(i) = yi 

X2 = Va(2) = y2 

= ya{3) = yi 
Xi = 2/0- (4) = ys 

X5 = ya{5) = 2/3- 



(9) 



Before we find the solutions we conclude that the number of solution is 

^nfc— 7 f^5-5 — 23 ^2 g 

First we have that 

0:4 = j/5 = 2 • 24 = 21 (mod 3^). 

We lift the result to Z243, that is, 

X4 = 2/5 = 21 + 3'*/i 21 + 3^/i, /leZd, =Z9. 

Hence, X4 G {21,48,75,102,129,156,183,210,237}. Since 36 is divisible by 9, 
neither 2/4 or 2/3 depends on 2/4, which follows from 

36 • (21 + 3^/1) = 36 • 21 EE 27 (mod 3^). 

We get that 2:5 = 2/4 = 179 and X3 = y^ = 32. The two first congruences in the 
system (|9]) gives us the nine solutions: 

(43,127,32,21,179) (151,73,32,48,179) (16,19,32,75,179) 
(124,208,32,102,179) (232,154,32,129,179) (97,100,32,156,179) 
(205,46,32,183,179) (70,235,32,210,179) (178,181,32,237,179). 
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6 Discussion 



In this paper we made the choice to you use as elementary methods as possible. 
Mainly because the problem is of an elementary nature but also because of the 
interest from applications. It is possible to lift both the problem and its solution 
into the context of free modules, see [S], over the ring Z/mZ. The derivation 
will when be almost identical but with different vocabulary. 

By small changes it is also possible to solve system of congruence equations 
with different modulus. We construct a system modulo the least common mul- 
tiple of all the moduli, that are equivalent to the original system. This is the 
same technique that is mentioned in ^y. 
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